Parentage analysis with genetic markers in natural populations. I. The expected proportion of offspring with unambiguous paternity.

Recent studies indicate that polymorphic genetic markers are potentially helpful in resolving genealogical relationships among individuals in a natural population. Genetic data provide opportunities for paternity exclusion when genotypic incompatibilities are observed among individuals, and the present investigation examines the resolving power of genetic markers in unambiguous positive determination of paternity. Under the assumption that the mother for each offspring in a population is unambiguously known, an analytical expression for the fraction of males excluded from paternity is derived for the case where males and females may be derived from two different gene pools. This theoretical formulation can also be used to predict the fraction of births for each of which all but one male can be excluded from paternity. We show that even when the average probability of exclusion approaches unity, a substantial fraction of births yield equivocal mother-father-offspring determinations. The number of loci needed to increase the frequency of unambiguous determinations to a high level is beyond the scope of current electrophoretic studies in most species. Applications of this theory to electrophoretic data on Chamaelirium luteum (L.) shows that in 2255 offspring derived from 273 males and 70 females, only 57 triplets could be unequivocally determined with eight polymorphic protein loci, even though the average combined exclusionary power of these loci was 73%. The distribution of potentially compatible male parents, based on multilocus genotypes, was reasonably well predicted from the allele frequency data available for these loci. We demonstrate that genetic paternity analysis in natural populations cannot be reliably based on exclusionary principles alone. In order to measure the reproductive contributions of individuals in natural populations, more elaborate likelihood principles must be deployed.

Altered transmission of HOX and apoptotic SNPs identify a potential common pathway for clubfoot.

Clubfoot is a common birth defect that affects 135,000 newborns each year worldwide. It is characterized by equinus deformity of one or both feet and hypoplastic calf muscles. Despite numerous study approaches, the cause(s) remains poorly understood although a multifactorial etiology is generally accepted. We considered the HOXA and HOXD gene clusters and insulin-like growth factor binding protein 3 (IGFBP3) as candidate genes because of their important roles in limb and muscle morphogenesis. Twenty SNPs from the HOXA and HOXD gene clusters and 12 SNPs in IGFBP3 were genotyped in a sample composed of non-Hispanic white and Hispanic multiplex and simplex families (discovery samples) and a second sample of non-Hispanic white simplex trios (validation sample). Four SNPs (rs6668, rs2428431, rs3801776, and rs3779456) in the HOXA cluster demonstrated altered transmission in the discovery sample, but only rs3801776, located in the HOXA basal promoter region, showed altered transmission in both the discovery and validation samples (P = 0.004 and 0.028). Interestingly, HOXA9 is expressed in muscle during development. An SNP in IGFBP3, rs13223993, also showed altered transmission (P = 0.003) in the discovery sample. Gene-gene interactions were identified between variants in HOXA, HOXD, and IGFBP3 and with previously associated SNPs in mitochondrial-mediated apoptotic genes. The most significant interactions were found between CASP3 SNPS and variants in HOXA, HOXD, and IGFBP3. These results suggest a biologic model for clubfoot in which perturbation of HOX and apoptotic genes together affect muscle and limb development, which may cause the downstream failure of limb rotation into a plantar grade position.

Detecting the signature of natural selection with microsatellites

Natural selection is one of the major factors in the evolution of all organisms. Detecting the signature of natural selection has been a central theme in evolutionary genetics. With the availability of microsatellite data, it is of interest to study how natural selection can be detected with microsatellites. ^ The overall aim of this research is to detect signatures of natural selection with data on genetic variation at microsatellite loci. The null hypothesis to be tested is the neutral mutation theory of molecular evolution, which states that different alleles at a locus have equivalent effects on fitness. Currently used tests of this hypothesis based on data on genetic polymorphism in natural populations presume that mutations at the loci follow the infinite allele/site models (IAM, ISM), in the sense that at each site at most only one mutation event is recorded, and each mutation leads to an allele not seen before in the population. Microsatellite loci, which are abundant in the genome, do not obey these mutation models, since the new alleles at such loci can be created either by contraction or expansion of tandem repeat sizes of core motifs. Since the current genome map is mainly composed of microsatellite loci and this class of loci is still most commonly studied in the context of human genome diversity, this research explores how the current test procedures for testing the neutral mutation hypothesis should be modified to take into account a generalized model of forward-backward stepwise mutations. In addition, recent literature also suggested that past demographic history of populations, presence of population substructure, and varying rates of mutations across loci all have confounding effects for detecting signatures of natural selection. ^ The effects of the stepwise mutation model and other confounding factors on detecting signature of natural selection are the main results of the research. ^

Antigenic variation in surface proteins of Trypanosoma cruzi (Chagas' disease)

Chagas' disease, a devastating illness in the Western Hemisphere, is caused by the protozoan parasite Trypanosoma cruzi. Transmission is via bloodsucking insect vectors, congenitally, or through blood transfusion and/or organ transplantation. A significant percentage of heart-related illnesses and deaths each year are attributable to the number of persons with Chagas' disease. Currently, there is no FDA-approved routine screening of the U.S. blood supply being conducted by blood banks. The only current commercial assays available for detection of Trypanosoma cruzi are based on South American isolates, which may differ antigenically from those found in the US. In this study, the assay used intact parasites as antigen in an ELISA-type assay. Therefore, serological differences presumably reflected variations in surface antigens. The basis of differential antibody binding to these antigens is unknown. In this study, biochemical characterization and genetic polymorphism analysis will be performed on three defined surface proteins of T. cruzi epimastigotes.^

An Sp1/Sp3 binding polymorphism confers methylation protection.

Hundreds of genes show aberrant DNA hypermethylation in cancer, yet little is known about the causes of this hypermethylation. We identified RIL as a frequent methylation target in cancer. In search for factors that influence RIL hypermethylation, we found a 12-bp polymorphic sequence around its transcription start site that creates a long allele. Pyrosequencing of homozygous tumors revealed a 2.1-fold higher methylation for the short alleles (P<0.001). Bisulfite sequencing of cancers heterozygous for RIL showed that the short alleles are 3.1-fold more methylated than the long (P<0.001). The comparison of expression levels between unmethylated long and short EBV-transformed cell lines showed no difference in expression in vivo. Electrophorectic mobility shift assay showed that the inserted region of the long allele binds Sp1 and Sp3 transcription factors, a binding that is absent in the short allele. Transient transfection of RIL allele-specific transgenes showed no effects of the additional Sp1 site on transcription early on. However, stable transfection of methylation-seeded constructs showed gradually decreasing transcription levels from the short allele with eventual spreading of de novo methylation. In contrast, the long allele showed stable levels of expression over time as measured by luciferase and approximately 2-3-fold lower levels of methylation by bisulfite sequencing (P<0.001), suggesting that the polymorphic Sp1 site protects against time-dependent silencing. Our finding demonstrates that, in some genes, hypermethylation in cancer is dictated by protein-DNA interactions at the promoters and provides a novel mechanism by which genetic polymorphisms can influence an epigenetic state.

Multiple independent genetic factors at NOS1AP modulate the QT interval in a multi-ethnic population.

Extremes of electrocardiographic QT interval are associated with increased risk for sudden cardiac death (SCD); thus, identification and characterization of genetic variants that modulate QT interval may elucidate the underlying etiology of SCD. Previous studies have revealed an association between a common genetic variant in NOS1AP and QT interval in populations of European ancestry, but this finding has not been extended to other ethnic populations. We sought to characterize the effects of NOS1AP genetic variants on QT interval in the multi-ethnic population-based Dallas Heart Study (DHS, n = 3,072). The SNP most strongly associated with QT interval in previous samples of European ancestry, rs16847548, was the most strongly associated in White (P = 0.005) and Black (P = 3.6 x 10(-5)) participants, with the same direction of effect in Hispanics (P = 0.17), and further showed a significant SNP x sex-interaction (P = 0.03). A second SNP, rs16856785, uncorrelated with rs16847548, was also associated with QT interval in Blacks (P = 0.01), with qualitatively similar results in Whites and Hispanics. In a previously genotyped cohort of 14,107 White individuals drawn from the combined Atherosclerotic Risk in Communities (ARIC) and Cardiovascular Health Study (CHS) cohorts, we validated both the second locus at rs16856785 (P = 7.63 x 10(-8)), as well as the sex-interaction with rs16847548 (P = 8.68 x 10(-6)). These data extend the association of genetic variants in NOS1AP with QT interval to a Black population, with similar trends, though not statistically significant at P<0.05, in Hispanics. In addition, we identify a strong sex-interaction and the presence of a second independent site within NOS1AP associated with the QT interval. These results highlight the consistent and complex role of NOS1AP genetic variants in modulating QT interval.

Population amalgamation and genetic variation: observations on artificially agglomerated tribal populations of Central and South America.

The interpretation of data on genetic variation with regard to the relative roles of different evolutionary factors that produce and maintain genetic variation depends critically on our assumptions concerning effective population size and the level of migration between neighboring populations. In humans, recent population growth and movements of specific ethnic groups across wide geographic areas mean that any theory based on assumptions of constant population size and absence of substructure is generally untenable. We examine the effects of population subdivision on the pattern of protein genetic variation in a total sample drawn from an artificial agglomerate of 12 tribal populations of Central and South America, analyzing the pooled sample as though it were a single population. Several striking findings emerge. (1) Mean heterozygosity is not sensitive to agglomeration, but the number of different alleles (allele count) is inflated, relative to neutral mutation/drift/equilibrium expectation. (2) The inflation is most serious for rare alleles, especially those which originally occurred as tribally restricted "private" polymorphisms. (3) The degree of inflation is an increasing function of both the number of populations encompassed by the sample and of the genetic divergence among them. (4) Treating an agglomerated population as though it were a panmictic unit of long standing can lead to serious biases in estimates of mutation rates, selection pressures, and effective population sizes. Current DNA studies indicate the presence of numerous genetic variants in human populations. The findings and conclusions of this paper are all fully applicable to the study of genetic variation at the DNA level as well.

CHARACTERIZING A NOVEL GENETIC LOCUS ASSOCIATED WITH FAMILIAL CO-OCCURRENCE OF THORACIC AORTIC ANEURYSMS AND INTRACRANIAL ANEURYSMS

The Mendelian inheritance of genetic mutations can lead to adult-onset cardiovascular disease. Several genetic loci have been mapped for the familial form of Thoracic Aortic Aneurysms (TAA), and many causal mutations have been identified for this disease. Intracranial Aneurysms (ICA) also show linkage heterogeneity, but no mutations have been identified causing familial ICA alone. Here, we characterized a large family (TAA288) with an autosomal dominant pattern of inherited aneurysms. It is intriguing that female patients predominantly present with ICA and male patients predominantly with TAA in this family. To identify a causal mutation in this family, a genome-wide linkage analysis was previously performed on nine members of this family using the 50k GenChips Hind array from Affymetrix. This analysis eventually identified a single disease-segregating locus, on chromosome 5p15. We build upon this previous analysis in this study, hypothesizing that a genetic mutation inherited in this locus leads to the sex-specific phenotype of TAA and ICA in this family First we refined the boundaries of the 5p15 disease linked locus down to the genomic coordinates 5p15: 3,424,465- 6,312,925 (GRCh37/hg19 Assembly). This locus was named the TAA288 critical interval. Next, we sequenced candidate genes within the TAA288 critical interval. The selection of genes was simplified by the relatively small number of well-characterized genetic elements within the region. Seeking novel or rare disease-segregating variants, we initially observed a single point alteration in the metalloproteinase gene ADAMTS16 fulfilling this criteria. This variant was later classified as a low-frequency population polymorphism (rs72647757), but we continued to explore the potential role of the ADAMTS16 as the cause of disease in TAA288. We observed that fibroblasts cultured from TAA288 patients consistently upregulated the expression of this gene more strongly compared to matched control fibroblasts when treated with the cytokine TGF-β1, though there was some variation in the exact nature of this expression. We also observed evidence that this protein is expressed at elevated levels in aortic aneurysm tissue from patients with mutations in the gene TGFBR2 and Marfan syndrome, shown by immunohistochemical detection of this protein.

Genetic variation in genes for the xenobiotic-metabolizing enzymes CYP1A1, EPHX1, GSTM1, GSTT1, and GSTP1 and susceptibility to colorectal cancer in Lynch syndrome.

Individuals with Lynch syndrome are predisposed to cancer due to an inherited DNA mismatch repair gene mutation. However, there is significant variability observed in disease expression likely due to the influence of other environmental, lifestyle, or genetic factors. Polymorphisms in genes encoding xenobiotic-metabolizing enzymes may modify cancer risk by influencing the metabolism and clearance of potential carcinogens from the body. In this retrospective analysis, we examined key candidate gene polymorphisms in CYP1A1, EPHX1, GSTT1, GSTM1, and GSTP1 as modifiers of age at onset of colorectal cancer among 257 individuals with Lynch syndrome. We found that subjects heterozygous for CYP1A1 I462V (c.1384A>G) developed colorectal cancer 4 years earlier than those with the homozygous wild-type genotype (median ages, 39 and 43 years, respectively; log-rank test P = 0.018). Furthermore, being heterozygous for the CYP1A1 polymorphisms, I462V and Msp1 (g.6235T>C), was associated with an increased risk for developing colorectal cancer [adjusted hazard ratio for AG relative to AA, 1.78; 95% confidence interval, 1.16-2.74; P = 0.008; hazard ratio for TC relative to TT, 1.53; 95% confidence interval, 1.06-2.22; P = 0.02]. Because homozygous variants for both CYP1A1 polymorphisms were rare, risk estimates were imprecise. None of the other gene polymorphisms examined were associated with an earlier onset age for colorectal cancer. Our results suggest that the I462V and Msp1 polymorphisms in CYP1A1 may be an additional susceptibility factor for disease expression in Lynch syndrome because they modify the age of colorectal cancer onset by up to 4 years.

Genetic networks controlling retinal injury.

PURPOSE: The present study defines genomic loci underlying coordinate changes in gene expression following retinal injury. METHODS: A group of acute phase genes expressed in diverse nervous system tissues was defined by combining microarray results from injury studies from rat retina, brain, and spinal cord. Genomic loci regulating the brain expression of acute phase genes were identified using a panel of BXD recombinant inbred (RI) mouse strains. Candidate upstream regulators within a locus were defined using single nucleotide polymorphism databases and promoter motif databases. RESULTS: The acute phase response of rat retina, brain, and spinal cord was dominated by transcription factors. Three genomic loci control transcript expression of acute phase genes in brains of BXD RI mouse strains. One locus was identified on chromosome 12 and was highly correlated with the expression of classic acute phase genes. Within the locus we identified the inhibitor of DNA binding 2 (Id2) as a candidate upstream regulator. Id2 was upregulated as an acute phase transcript in injury models of rat retina, brain, and spinal cord. CONCLUSIONS: We defined a group of transcriptional changes associated with the retinal acute injury response. Using genetic linkage analysis of natural transcript variation, we identified regulatory loci and candidate regulators that control transcript levels of acute phase genes.

HLA-DPB1 and DPB2 are genetic loci for systemic sclerosis: a genome-wide association study in Koreans with replication in North Americans.

OBJECTIVE: To identify systemic sclerosis (SSc) susceptibility loci via a genome-wide association study. METHODS: A genome-wide association study was performed in 137 patients with SSc and 564 controls from Korea using the Affymetrix Human SNP Array 5.0. After fine-mapping studies, the results were replicated in 1,107 SSc patients and 2,747 controls from a US Caucasian population. RESULTS: The single-nucleotide polymorphisms (SNPs) (rs3128930, rs7763822, rs7764491, rs3117230, and rs3128965) of HLA-DPB1 and DPB2 on chromosome 6 formed a distinctive peak with log P values for association with SSc susceptibility (P=8.16x10(-13)). Subtyping analysis of HLA-DPB1 showed that DPB1*1301 (P=7.61x10(-8)) and DPB1*0901 (P=2.55x10(-5)) were the subtypes most susceptible to SSc in Korean subjects. In US Caucasians, 2 pairs of SNPs, rs7763822/rs7764491 and rs3117230/rs3128965, showed strong association with SSc patients who had either circulating anti-DNA topoisomerase I (P=7.58x10(-17)/4.84x10(-16)) or anticentromere autoantibodies (P=1.12x10(-3)/3.2x10(-5)), respectively. CONCLUSION: The results of our genome-wide association study in Korean subjects indicate that the region of HLA-DPB1 and DPB2 contains the loci most susceptible to SSc in a Korean population. The confirmatory studies in US Caucasians indicate that specific SNPs of HLA-DPB1 and/or DPB2 are strongly associated with US Caucasian patients with SSc who are positive for anti-DNA topoisomerase I or anticentromere autoantibodies.

The angiotensin-converting enzyme insertion/deletion polymorphism and its associations with carotid artery wall thickness and coronary heart disease: The atherosclerosis risk in communities study

Coronary heart disease (CHD) is the leading cause of death in the United States. Recently, renin-angiotensin system (RAS) was found associated with atherosclerosis formation, with angiotensin II inducing vascular smooth muscle cell growth and migration, platelet activation and aggregation, and stimulation of plasminogen activator inhibitor-1. Angiotensin II is converted from angiotensin I by angiotensin I-converting enzyme (ACE) and this enzyme is mainly genetically determined. The ACE gene has been assigned to chromosome 17q23 and an insertion/deletion (I/D)polymorphism has been characterized by the presence/absence of a 287 bp fragment in intron 16 of the gene. The two alleles form three genotypes, namely, DD, ID and II and the DD genotype has been linked to higher plasma ACE levels and cell ACE activity.^ In this study, the association between the ACE I/D polymorphism and carotid artery wall thickness measured by B-mode ultrasound was investigated in a biracial sample, and the association between the gene and incident CHD was investigated in whites and if the gene-CHD association in whites, if any, was due to the gene effect on atherosclerosis. The study participants are from the prospective Atherosclerosis Risk in Communities (ARIC) Study, including adults aged 45 to 65 years. The present dissertation used a matched case-control design for studying the associations of the ACE gene with carotid artery atherosclerosis and an unmatched case-control design for the association of the gene with CHD. A significant recessive effect of the D allele on carotid artery thickness was found in blacks (OR = 3.06, 95% C.I: 1.11-8.47, DD vs. ID and II) adjusting for age, gender, cigarette smoking, LDL-cholesterol and diabetes. No similar associations were found in whites. The ACE I/D polymorphism is significantly associated with coronary heart disease in whites, and while stratifying data by carotid artery wall thickness, the significant associations were only observed in thin-walled subgroups. Assuming a recessive effect of the D allele, odds ratio was 2.84 (95% C.I:1.17-6.90, DD vs. ID and II) and it was 2.30 (95% C.I:1.22-4.35, DD vs. ID vs. II) assuming a codominant effect of the D allele. No significant associations were observed while comparing thick-walled CHD cases with thin-walled controls. Following conclusions could be drawn: (1) The ACE I/D polymorphism is unlikely to confer appreciable increase in the risk of carotid atherosclerosis in US whites, but may increases the risk of carotid atherosclerosis in blacks. (2) ACE I/D polymorphism is a genetic risk factor for incident CHD in US whites and this effect is separate from the chronic process of atherosclerosis development. Finally, the associations observed here are not causal, since the I/D polymorphism is in an intron, where no ACE proteins are encoded. ^

Global patterns of DNA polymorphism in noncoding regions in human populations

DNA sequence variation is currently a major source of data for studying human origins, evolution, and demographic history, and for detecting linkage association of complex diseases. In this dissertation, I investigated DNA variation in worldwide populations from two ∼10 kb autosomal regions on 22q11.2 (noncoding) and 1q24 (introns). A total of 75 variant sites were found among 128 human sequences in the 22q11.2 region, yielding an estimate of 0.088% for nucleotide diversity (π), and a total of 52 variant sites were found among 122 human sequences in the 1q24 region with an estimated π value of 0.057%. The data from these two regions and a 10 kb noncoding region on Xq13.3 all show a strong excess of low-frequency variants in comparison to that expected from an equilibrium population, indicating a relatively recent population expansion. The effective population sizes estimated from the three regions were 11,000, 12,700, and 8,600, respectively, which are close to the commonly used value of 10,000. In each of the two autosomal regions, the age of the most recent common ancestor (MRCA) was estimated to be older than 1 million years among all the sequences and ∼600,000 years among non-African sequences, providing first evidence from autosomal noncoding or intronic regions for a genetic history of humans much more ancient than the emergence of modern humans. The ancient genetic history of humans indicates no severe bottleneck during the evolution of humans in the last half million years; otherwise, much of the ancient genetic history would have been lost during a severe bottleneck. This study strongly suggests that both the “out of Africa” and the multiregional models are too simple for explaining the evolution of modern humans. A compilation of genome-wide data revealed that nucleotide diversity is highest in autosomal regions, intermediate in X-linked regions, and lowest in Y-linked regions. The data suggest the existence of background selection or selective sweep on Y-linked loci. In general, the nucleotide diversity in humans is low compared to that in chimpanzee and Drosophila populations. ^

Are functional and genetic components of platelets linked to coronary artery disease: A case-control study

Coronary artery disease (CAD) is a multifactorial disease process involving behavioral, inflammatory, clinical, thrombotic, and genetic components. Previous epidemiologic studies focused on identifying behavioral and demographic risk factors of CAD, but none focused on platelets. Current platelet literature lacks the known effects of platelet function and platelet receptor polymorphisms on CAD. This case-control analysis addressed these issues by analyzing data collected for a previous study. Cases were individuals who had undergone CABG and thus had been diagnosed with CAD, while the controls were volunteers presumed to be CAD free. The platelet function variables analyzed included fibrinogen Von Willebrand Factor activity (VWF), shear-induced platelet aggregation (SIPA), sCD40L, and mean platelet volume; and the platelet polymorphisms studied included PIA, α2 807, Ko, Kozak, and VNTR. Univariate analysis found fibrinogen, VWF, SIPA, and PIA to be independent risk factors of CAD. Logistic regression was used to build a predictive model for CAD using the platelet function and platelet polymorphism data adjusted for age, sex, race, and current smoking status. A model containing only platelet polymorphisms and their respective receptor densities, found polymorphisms within GPIbα to be associated with CAD, yielding an 86% (95% C.I. 0.97–3.55) increased risk with the presence of at least 1 polymorphism in Ko, Kozak, or VNTR. Another model included both platelet function and platelet polymorphism data. Fibrinogen, the receptor density of GPIbα, and the polymorphism in GPIa-IIa (α2 807) were all associated with CAD with odds ratios of 1.10, 1.04, and 2.30 for fibrinogen (10mg/dl increase), GPIbα receptors (1 MFI increase), and GPIa-IIa, respectively. In addition, risk estimates and 99% confidence intervals adjusted for race were calculated to determine if the presence of a platelet receptor polymorphism was associated with CAD. The results were as follows: PIA (1.64, 0.74–3.65); α2 807 (1.35, 0.77–2.37); Ko (1.71, 0.70–4.16); Kozak (1.17, 0.54–2.52); and VNTR (1.24, 0.52–2.91). Although not statistically significant, all platelet polymorphisms were associated with an increased risk for CAD. These exploratory findings indicate that platelets do appear to have a role in atherosclerosis and that anti-platelet drugs targeting GPI-IIa and GPIbα may be better treatment candidates for individuals with CAD. ^